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Sfcbstract 

S^Jfconsistent equations which describe the order parameter and chemical potential be- 
haviour in 2D metalic system as functions of external magnetic field, H, temperature, T, 
anH carrier density, n, are obtained. It is shown that for the case of the local pairs (low 
^jrier density and negative chemical potential fi) the derivative dH C2 (T)/dT at T = T c is 
essentially less then that for the system with Cooper pairs (high carrier density and positive 
comical potential u. ^> T c ). It is found that in magnetic field satisfying the quantum limit 
psjjterium the system is characterized by non-trivial inhomogeneuos order parameter which 
Sh exist at rather high temperatures. 

\o : 

8 Introduction 

I^spite of the great deal of efforts directed toward solving the problem, the nature of high 
t^nperature superconductivity (HTSC) retains to be unknown. However, those efforts have 

! been done in vain. At present one can observe a common consent about the most 
cteiracteristic properties of HTSC [I], |2[. Such properties are: a) quasi-two-dimensional 
tQiaviour of the conductivity in the normal phase; b) relatively low (at least considerably 
lo^er than in ordinary metals) density of carriers. There are also many other distinctive 
properties which may appear to be very important for understanding the nature of HTSC. 
!$Ht! it is unlikely to describe all the characteristic features of HTSC in the framework of the 
oMy theoretical approach. So, very often the role of one of them is studied and then the 
total picture is created by means of summing up all the results in some "artistic" way. 

One of the very important and intensively studied questions in the theory of supercon- 
ductivity (SC) is the question about a dependence of SC properties of a system on the carrier 
density. The first who emphasized on the relevance of such a problem was Legget ||. Later 
this problem was also discussed in Refs. [|], ||, f| . In Ref. @ the existence of a crossover driven 
by the carrier density was clearly established. In that paper the authors have considered 
a three dimensional model with a local attractive interaction and shown that the gradual 
transition from low to high carrier density is accompanied by a crossover from local pair SC 
to the SC with Cooper pairs. As for HTSC, an intermidiate case seems to be realized in it. 

In this paper we continue the analysis of the papers 0] and H. 
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2 Model and general discussion 



As in the papers H §, we shall focus our attention only on the qualitative side of the 
problem. So, for our purposes, we can take the simplest Hamiltonian which describes a 
system of charge carriers (fermions) with a local attraction between them: 



H 
H(t) 



d 2 rH(v 



— ^— (dj — i-Aj 
2m V 3 c 3 



9. 



(2) 



where Aj is a vector potential corresponding to a magnetic field, H, (note that in two 
dimensions magnetic field is pseudoscalar function), e and m are charge and mass of the 
carriers, respectively, g > is an attractive coupling between these carriers and /i is the 
chemical potential. 



Introducing Nambu notations for fermion field ^ 
more convenient form: 



we can rewrite 



W(r) 



V 2 



\ 2m / 



m a 



(3) 



where r 3 , r + = (n + zr 2 )/2, r_ = (rj — ir 2 )/2 are Pauli matrices and the covariant derivative 
has the form X> 5 = 9j — ier 3 Aj/c. Then the partition function is expressed through the 
Hamiltonian as: 

Z = exp(-Sl/T) = Trexp(-HfT), (4) 

Since the thermodynamical potential, Q = Q(V,T, /x), is the function of the chemical po- 
tential (besides the volume and the temperature) and we are interested in a dependence of 
all physical values on the density of the carriers, it is necessary to write down the second 
equation which links the density, n, and the chemical potential, \i: 

(5) 



n 



Vdfi' 



The expression for the partition function with the Hamiltonian (^) can be represented in a 
path integral form: 



Z = J [d^d$>] exp - / dr j d 2 r (V<9 T ^ + H{r] 



1 

T' 



(6) 



where Grasmann variables \l/(r;r) satisfy the antiperiodic boundary condition \l/(0;r) = 
—^>(j3;r). After usual introduction of an auxiliary scalar field by means of the Habbard- 
Stratanovich trick: 



exp — J dr J d 



+ ^ 



d T - r 3 



^2 

2m 



l$l 



+ 



(7) 
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we can perform (at least formally) the integration over Grasmann variables and represent 
the result through the effective action: 

S ef f ($,$*) = -TrLnG- 1 + - J J d 2 r|$| 2 , (8) 

3 o 

depending only on the scalar field $. However, the next integration over scalar field we can 
perform only approximately, for example, using so called saddle point formalism. In our 
problem the "saddle point" is defined by the equation: 

Z\r-r) =MG(T,r;r,r)r + ] + - = 0, (9) 

where Green's function, G, of interacting fermions is defined as the solution of the equation: 
V 2 



G(n, r 2 ; ri, r 2 ) = <y(n - r 2 )5(ri - r 2 ), (10) 



with the boundary condition G(t\ + (3, r 2 ; ri, r 2 ) = —G(ti, r 2 ; ri, r 2 ). 

Let $ be a solution to the equation @, then the partition function in the next to the 
leading approximation takes a form of the Gauss type path integral: 

Z = exp(-S e ff) J[d®d®*]exp - J dn J dr 2 



■ J d 2 r, J d 2 r 2 $*(r 1 ;r 1 )r- 1 (ri,r 2 ;r 1 ,r 2 )$(r 2 ;r 2 )] > (11) 
where S e ff = 5 e //( < l>, $*). New field $ describes fluctuations and T is its propagator: 
r _1 (ri,r 2 ;ri,r 2 ) = -5(n - r 2 )S(r 1 - r 2 ) + 

+ tr G(ri,r 2 ;ri,r 2 )r_G(r 2 ,r 1 ;r2,r 1 )r + . (12) 
L J $ 

For the partition function in this approximation we have 

Z = exp(-S eff - TrLnY- 1 ). (13) 

So, as follows from Eqs.(|4]) and (|5|), the expression for the carrier density takes the form of 
a sum of two different terms: 

n = ^Tr[r 3 G] - ^^{TrLnY~ l ). (14) 

The first term is expressed through the fermion propagator and the second through the 
propagator of scalar fluctuations. Having such a representation, we shall refer to the first 
term as fermion part in density of the carriers and to the second as boson one. The fact that 
the ratio of fermions and composite bosons (fluctuations can also be interpreted as a field of 
two-fermion composite particles) is determined by the dynamics is the main virtue of this 
model. 

At this point we have, in principle, a closed selfconsistent system of two equations which 
completely describe a dependence of the order parameter and the chemical potential as 
functions of the temperature and carrier density. 
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3 On the critical line equation 



As was indicated in the previous section Eqs.@ and ([L4]) completely describe the behaviour 
of the order parameter and the chemical potential of the system as a function of "external" 
parameters T, B and n. But analizing these equations in general case is an unsolvable 
problem. So we restrict ourselves only to analizing the behaviour of the system near the 
critical line (related to the second type phase transition). Such a choice simplifies the problem 
considerablly since when the system is in the nearcritical region we have a natural small value 
(order parameter) and as a consequence we can apply pertubation theory in this value. 

In order to simplify the problem further we assume that the solution to the equation 
(§) does not depend on "time" coordinate r. In this follows from Eq.(|i~0|), Green's 

function depends only on the difference of "time" variables (t% — r 2 ). So, after taking into 
account the boundary conditions, Green's function can be expanded into the Fourier series: 

oo 

G(n - r 2 ; r x , r 2 ) = T G n{*i, r 2 ) exp[-iu n (n - r 2 )], w„ = ixT(2n + 1). (15) 

n=— oo 

Solving the equation for the Green's function (fl~0| ) in linear approximation in order parameter 
and substituting it into (H), we come to the following integral equation: 

$(r) = J rf 2 r 1 K 1 (r,r 1 )$(r 1 ), (16) 



K 1 (r,r 1 ) = -gT £ tr Gf(r, n)r_Gf( ri , r)r. 



(17) 



where Gffl(r, ri) denotes the Green's function at $ = which can be easily found using the 
Schwinger proper time method ||10|| . Here we write down this function without deriving (for 
details see appendix A in Ref.||ll||): 



Gftn.rJ = e W (-^T 3 (x iy2 - yi x 2 )jG^"Hyi-r 2 ), (18) 



Strictly speaking, Eq.flTTD which has the form of a spectral problem is valid only on the 
critical line where the order parameter equals zero. The condition indicating the existence 
of a solution to this equation determines the critical line in the problem. If we were going to 
leave the critical line, we would need to take into account the next nonlinear term. In this 
paper we are interested only in the behaviour on the critical line. 

After substituting the explicit expression for the Green's function into (|T7|) we get: 

$(r) = Jd 2 r 1 Kl hom) (r-r 1 )exp^(xy l -yx l )^(r 1 ), (20) 
where / is the magnetic length and the homogeneous kernel, x[ hom \ determined by 



i^ om V,ri) = -gT £ tr G^ om \v)T„G^ m \*)T + 



(21) 
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As was shown in |I2" , |13f , Eq. fl2"0|) has the exact solution: 

r.2 \ 



A exp 



r 

2Z2 



where in the simplest case A is constant. After substituting 
condition determining the critical line: 



(22) 

into (^U|) we come to the 



gTf 
(2vr) ; 



c>o r r 

Y, d 2 Kd 2 kexp 



2l2- 



tr 



^(hom) 



K 



K 



(23) 



The second equation which links the density and the chemical potential is the following: 



T d 



n 



(27T) 1 



£ / d 2 ktr[r 3 G^ m \k)] - ^(Trkr - :l ). 



(24) 



where To is defined by the expression similar to ([121) but with Green's functions at $ = 0. 

Eqs.(|23"D and ([24]) selfconsistently describe the critical parameters u(T) and H C2 (T) at 
every given carrier density. Below we shall consider two limiting cases of low and high 
magnetic fields. 



4 Low magnetic field 

This case is realized when the cyclotron frequency is much lower then the temperature, 
u H < T c . Expanding Eqs.flU and flU) in power seriers in I 1 and keeping only the lowest 
term in magnetic field, we get the equation: 



Air 
gm 



w 



du u — u 4T x , 
tanh — — 2^ 



x U '"'du d u 2 (u + fi) 



u — jl 



2T 



mf 



n=0 



u du ioi + u 



2 ' 



(25) 



where W is the band width. 

The critical temperature, T c , when the field is swiched off is defined by the equations: 



gm 



du u- u c 

tanh ■ 



u - u c 



2T r 



/i c = /i(T c ), 



(26) 



which have been analyzed in detail in Ref. ||. Substituting the expression for the coupling 
(p6|) into (|25"D and taking the limit W — > oo, we come to: 



du 



tanh(w — u)/2T tanh(w — u c )/2T c 



u — u 



u— u c 



2T ~ 

777/2 



1 ttu \x arctan(/i/u; ri ) 



2ujI 



(27) 
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Approaching the critical line at low but nonzero magnetic field, we should come to the value 
of the temperature close to T c . So, expanding the left hand side of the Eq. fl27D in (T — T c ) 
we can find the slope of the critical line at T c : 

e ( dH C2 \ 1 + tanh(/i c /2T c ) - (fi c /T c )(d/i c /dT c ) tanh(/i c /2T c ) 

■ 



V /Tc l + [7£(3)/2ir*](jJL c /T c ) + (\(jL c \/2T e ) J {du/v?){u - tanhw) 

o 

Even though the slope of the critical line is given by the explicit expression fl28|) this does not 
solve the problem completely because, as was indicated above, we should add also the second 



equation (24) linking the chemical potential with the carrier density. From that equation 
we should also find the derivative <9/i c /<9T c = (<9///<9T)r c which is calculated at constant n. 
Since the analysis of the second equation for arbitrary fi is quite a difficult problem we shall 
consider only limiting cases. 

1. Local pairs. At first we consider the case of low density of carriers and strong interaction 
between fermions, i.e. ep <C |e&| where |e&| = 2W / exp(— An/rag) is the energy of bound state 
in two particle problem and ep = rnr/m. Analyzing the propagator of the scalar field 
describing fluctuations we come to the conclusion that in the lowest approximation in Z -1 
the propagator should be taken at zero magnetic field (for details see WW). Assuming that 



/j, c < and |/x c | ^> T c (the equivalence of these conditions and the inequality stated above 
can be easily established), it is not difficult to find the expression for carrier density in this 
limiting case: 



2mT c 
n = In 

7T 



exp 



2/i c - e b 



T, 



c 




O exp -!£! , (29) 



where we omitted exponentially small terms. 

In principle, Eqs.fl2"Bp and fl2"9p completely describe the critical line slope at T c , if the 
critical temperature at zero field is nonzero and the density of carriers is a finite constant. 
However, it can be easily shown that equations fl2"B|) and ( ^9]) are consistent only at T c = 0. 
This fact is in full agreement with the general statements about the role of fluctuations in 
low dimentional models with short-range interaction where, as is known, there is no room 



for a nontrivial order parameter |fL4|| . This problematic situation is not actual in real HTSC 
since they are not strictly two-dimensional, they are only quasi-two-dimentional. So, we 
come to the place where we are not able to use a strictly two-dimensional model. In order 
to avoid this obstacle we assume that the model will be treated correctly if we take into 
account quasi-two-dimensional character of the model only into the equation for carrier 
density which is extremely sensitive to the number of dimensions and leave the equation 
(P8|) without changes. 

Since in three dimensions the density is measured in different units, our enlargement 
of the phase volume of the system should be accompanied by a redefinition of the "two- 
dimentional" density in the form: n 2D = v? D K~ x where n 3D is the density found in the 3D 
anisitropic model and K z is some characteristic momentum. Chosing K^/4M ~ T where M 
is an effective mass in z-direction, we get: 

2mT 7 , u 1 ' 2 , . 

n = — ^— / du F ; — — -. (30) 

7t 2 J exp[w + (e - 2/i)/T - 1] K ' 
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Now we find the slope of the critical line without problems: 

e / dH C2 \ ( fi c \ 2 ( \/i c \\ , , \e b 







2T C r 


| 



4mc\dT J- \TJ — V 2TJ' ™ ~T Tc ~ ^ (31) 
Obtained expression shows that the slope of the critical line in the case of composite bosons 



is very small what is in qualitative agreement with the result of the paper |T5[ 



2. Cooper pairs. When the density of carriers is high or the attraction between fermions 
is relatively weak, i.e. ep 3> C6> we come to the SC of Cooper pairs. In terms of chemical 
potential this condition is equivalent to the inequality ji c ^> T c . Since the fluctuations are 
suppressed in any three-dimensional model and the mean field approximation is reliable 
when the density of carriers is high enough, we can assume that the same statement is 
correct in our "quasi-two-dimensional" model (though it cannot be shown in a direct way). 
With such an assumption the bosonic term in Eq . Q24|) can be omitted and we get: 



^^(l + exp^)-^, (32) 



and for the slope of the critical line: 

e (dH C2 \ 2tt 2 T c _ 7 



yUc ~e F >T c = -j2\e b \e F , (33) 

7T v 



where 7 is the Euler constant. 

So, the results received in this section show that the behaviour of the critical line, in 
particular its slope at T c essentially depends on the density of carriers in the system. At 
low density the slope is small and with increasing the density the slope becomes bigger (all 
other parameters are kept constant). Note that after achieving some maximum value the 
slope turns to decreasing since in the regime of Cooper pairs it decreases with increasing of 
the carrier density. 



5 Strong magnetic field 

When a magnetic field is so strong that it satisfies the condition, u>h 3> T c , we come to the 
so called quantum limit. It is obvious that in a very strong magnetic field, only the lowest 
Landau level in fermion spectrum plays an important role. So, omitting all those terms in 
fermion Green's function which correspond to the higher Landau levels, we get: 

r(hom) (]] _ 2exp(-k 2 / 2 ) 

This, in its turn, means that the density of states at the lowest level should be higher than 
the density of carriers in the system. It is also clear that the higher density of carriers is in 
the system, the stronger magnetic field should be applied in order to achieve the quamtum 
limit. The second condition when the approximation of the lowest Landau level is reliable 
is connected with the value of coupling constant. When the attraction between fermions is 
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very strong the substitution of the approximate expression ([53) instead of an exact Green's 
function will not be satisfactory since any strong interaction could considerably change the 
spectrum. 

At last, all analysis in this section is performed in the framework of the mean field ap- 
proximation. We beleive that the consideration taking into account fluctuations in quantum 
limit would have been worth-while problem, but as is easily seen from all our consideration, 
at first one need to find some simple formalism for treating quasi-two-dimensional models 
which is free from the obstacles incorporated into all two-dimensional models. It could be 
more realistic to consider an anisotropic 3D model, however as we can judge such a problem 
is much more difficult than ours and could be solved only numerically. 

Since the expression for Green's function is so simple, we come at once to the equation 
describing the critical line: 

rp _ - 2e F 

2h\{u H /\e b \)h\{u H /e F - 1) ' 

where we assumed that n < eH/nc (or e F < ojh). The obtained result is similar to the result 
of paper |16||. 



6 Conclusion 

In this paper we derived the system of selfconsistent equations describing second type phase 
transition in the simplest "quasi-two-dimensional" system with local attraction between 
fermions in an external magnetic field. We have shown that when the density of carriers 
in the system is relatively low the equation, linking the chemical potential and the carrier 
density, is nontrivial (// ^ e F ) and, as a result, this is reflected in the behaviour of the critical 
line (its slope at T c is small). In the opposite case of figh density of carriers the system can 
be described in the framework of the mean field approximation (as in BCS model). 

Our results are in qualitative agreement with the results obtained in [15j for composite 
Bose particles though analytical behaviour is different. 

The research was supported in part by the grant No.K5O100 of the International Science 
Foundation and by the International Soros Science Education Program (ISSEP) through the 
grant No.PSU052143. 
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